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Abstract: 

Top quark-antiquark (tt) pairs will be produced copiously at the Tevatron collider and in 
huge numbers at the LHC. This will make possible detailed investigations of the properties 
and interactions of this quark flavor. The analysis and interpretation of future data requires 
precise predictions of the hadronic production of tt pairs and of their subsequent decays. 
In this talk the reactions pp.pp — ?f + X — > +X are considered and results are 
presented of our calculation[|T]] of the dilepton angular distribution at next-to-leading order 
QCD, keeping the full dependence on the spins of the intermediate tt state. The angular 
distribution is determined for different choices of reference axes that can be identified 
with the t and t spin axes. While the QCD corrections to the leading-order distribution 
turn out to be small in the case of the LHC, we find them to be sizeable in the case of the 
Tevatron and find, moreover, the angular distribution to be sensitive to the parton content 
of the proton. 
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Speaker at the conference 



So far the top quark, which was discovered six years ago[|^] is still a relatively unex- 
plored particle, as compared with other quark flavors. This will change once the upgraded 
Tevatron collider and, in several years, the LHC will be in full operation. It is expected 
that about 10"* top quark- antiquark (tt) pairs per year will be produced at the Tevatron and 
more than 10^ tt pairs per year at the LHC. These large data sample will make feasible 
precise investigations of the properties and interactions of top quarks. 

An important aspect of top quark physics will then come into play, namely, spin 
physics with top quarks. It is well-known by now that the t quark is unique among the 
quark flavors in that it allows to study effects associated with its spin in a direct and unam- 
biguous way. This is due to its extremely short lifetime that prevents the top quark from 
forming hadronic bound states. It behaves like a highly instable "bare" quark. Therefore 
phenomena associated with the spins of the top quark and antiquark are reflected directly 
in the distributions and in the corresponding angular correlations of the jets, W bosons, or 
leptons into which the t and t decay. These distributions and correlations reflect the t and 
i polarizations and spin correlations which in turn characterize the t and t production and 
decay mechanism(s). 

On the theoretical side the spin correlations of hadronically produced tt pairs were 
studied some time ago[^ ^, |^ to leading order in the coupling a? of Quantum Chro- 
modynamics (QCD). There exists also an extensive literature, for example refs.[Q] and 
references therein, on how to exploit top-quark spin phenomena at hadron colliders in 
the search for new interactions. For instance the spin and parity of a new heavy reso- 
nance that strongly couples to tt in the s channel could be pinned down with spin correla- 
tions. With appropriate observables that reflect the t and/or i spins one can, for instance, 
check whether or not the V — A law is valid also for top decay, t Wb, or search for 
non-standard CP violation in ti production and/or decay. A prerequisite of this kind of 
experimental analysis is that these spin effects must be known as precisely as possible 
within the standard model (SM) of particle physics. Therefore we have determined the 
production of ti pairs by qq annihilation, gluon-gluon fusion and gluon-^(^) scattering at 
next-to-leading order (NLO) in the QCD coupling for arbitrary t and i spins [^. More- 
over, we have analyzed the hadronic production of pairs and their subsequent decays to 
order a^, keeping the full information on the spins of the intermediate tt state. Within the 
SM, where the main top-quark decay modes are t — > bW bqq'.bivg, the most powerful 
analyzers of the polarization of the top quark are the charged leptons, or the jets that orig- 
inate from quarks of weak isospin —1/2 produced by the decay of the W boson. Here we 
restrict ourselves to the channels where both t and t decay semileptonically, 

pp,pp^it+X ^i+£'-+X, (1) 

(i = e,^,x), and we present our predictions of the dileptonic angular distribution[jl|] that 
encodes the ti spin correlations. 

At the parton level the NLO analysis involves the following subprocesses: 

gg,qq bbi+i'-veve, (2) 

gg,qq ^ bbi+i'-veve>+g, (3) 
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g + q{q) ^ bM+i'-VeVe+q{q). 



(4) 



At the Tevatron tt production is dominated by quark-antiquark annihilation while at the 
LHC it is mainly due to gluon-gluon fusion. 

Because the total width Ft of the top quark is much smaller than its mass nit (^t /mt = 
0(1%)), we can expand the amplitudes of the parton reactions (|^) - around the poles 
of the unstable t and f quarks and keep only the leading term of this expansion, i.e., 
the residue of the double poles. The radiative corrections to the respective lowest-order 
amplitudes can be classified into so-called factorizable and non-factorizable corrections. 
We take into account the factorizable corrections to the above reactions for which the 
squared matrix element is of the form (''^^p °^ Tr[p(-^)i?("'')p(-^ )]. Here X = 1, ...,6 
labels the 6 amplitudes of (Q) - (Q), and r("-') denotes the respective spin density matrix 
for the production of on-shell ?f pairs. The superscript a labels the initial state and i lables 
the intermediate state, i.e., i = tt,ttg,ttq,ttq state. The decay density matrix p(-^)(p(^)) 
describes the normalized angular distribution of the decay of a polarized t{t) quark into 
+ anything in the rest frame of the t{t) quark. Note that for the reactions (|^) the 
squared matrix elements | fAf ^'^^ p have, for each X, three different contributions of the form 
because the final-state gluon is either associated with the tt production, 
the t, or t decay amplitude. 

Let us first discuss the production density matrices R^"''^ and the QCD-induced spin 
correlation effects at the level of the tt states. The spin density matrix is defined in 
terms of the respective transition matrix element T {a ^ i) . For instance, for gg tt we 
have 

4a'',pp' = T^ E m^\gg){gg\^'\t^h'). (5) 

^gg colors 
initial spins 

where the factor Ngg averages over the spins and colours of the initial pair of partons. The 
matrix structure of the R^'^'^'^ is (for ease of notation we drop the superscripts) 



+ Q/(o^)aa'((5')pp' , (6) 

The function A = Tr(7?)/4 determines the differential cross section with tt spins summed 
over. Because of parity invariance the vectors B, B can have, within QCD, only a compo- 
nent normal to the scattering plane. This component, which amounts to a normal polar- 
ization of the t and i quarks, is induced by the absorptive part of the scattering amplitude 
which, for the i = tt intermediate state, is non-zero to order . The normal polarization 
is quite small, both for ti production at the Tevatron and at the LHC [§]. The functions 
Cki encode the correlation between the t and t spins. 

In the computation of the R^"''^ to order we used dimensional regularization to 
treat both the ultraviolet and the infrared/coUinear singularities. Renormalization was 
performed using the MS prescription for the QCD coupling a? and the on-shell definition 
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of the top mass nit. For a = qq.gg the soft and coUinear singularities in R^'''"\ which 
appear as single and double poles in £ = (4 — D)/2, are cancelled after including the con- 
tributions from R^"'"s^ in the soft and coUinear limits and after mass factorization. For the 
latter we used the MS factorization scheme. The order a, production density matrices of 
(0) contain initial state coUinear singularities which are also removed by mass factoriza- 
tion. The soft and coUinear singularities were extracted by employing a simplified version 
of the phase-space slicing technique. 

From these density matrices one can obtain, in particular, the total parton cross sec- 
tions 6a(f) for the reactions gg,qq, q/{q)g ^ tt + X at NLO. We have computed these 
cross sectionslQ] and found excellent agreement with previous results[^ |T0|]. 

In order to study the tt spin correlations at the parton level we consider the following 
class of observables 

0=4(a-s,)(b-Sf) (7) 

where Sf,Sf denote the t and f spin operators, and a and b are reference directions that 
serve as spin axes. We choose the following vectors: 

a = , b = kf (helicity basis) , 
a = p, b = p (beam basis), 

a = df , b = df (off-diagonal basis). (8) 

Here k;(kf) denotes the direction of flight of the t{t) quark in the parton center-of-mass 
frame, and p is the unit vector along one of the hadronic beams in the laboratory frame. 
Furthermore d? is the axis with respect to which the spins of t and t produced by qq 
annihilation are 100 % correlated[^ [Tl[] to leading order in a^. (For gg tt one can 
show that no spin basis with this property exists.) 



Figure 1: The unnormalized expectation value of the spin correlation observable (0) in 

the beam basis. The plots show g^a^ (rj) (dotted line), gcP (ri)(full line) and g^}^ (rj) (dashed 
line) as a function of r\ for the qq (a), and gg (b) initial state. 
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The expectation values of the above observables O signify the degree of correlations 
among the ti spins when using the above spin quantization axes. If one identifies the MS 
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renormalization scale jjR with the mass factorization scale jup, /Jr= /Jf = fJ, and neglects 
all quark masses except for , then one can express the unnormalized expectation value 
of a spin-correlation observable in terms of dimensionless scaling functions as follows: 
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m. 



where a = gg, qq,q/ {q)g, s is the parton center-of-mass energy squared, and T] = ^ — 1 . 

As an example the functions ^i^^(Ti), g'^\v[) and gcP{T\) are shown in Fig|l| for (a = 
qq^gg) and the observable in the beam basis. The contributions of the q/ {q)g initial states 
to the spin correlations at the hadron level are very small. The results for the other spin 
observables are given in ref.[^. 

The QCD corrections to the unnormalized spin correlations are large close to thresh- 
old. This behaviour is due to the factor Oa which, in this order of perturbation theory, is 
non-zero at threshold due to Coulomb attraction. For hadronic observables these effects 
are damped by the parton distribution functions and integration over the momenta of the 
partons in the initial state. 

Next we discuss t and t decay. Here only the semileptonic decays of a polarized t{t) 
into ^^{^'^) -\- anything are considered. The decay density matrix p(-^)(p(-^')) has the form 
•^Pa'cc ~ (11 + i^+CJ ■ Q+)a'a where q+ describes the direction of flight of in the rest 
frame of the t quark and a, denote the Pauli matrices. The decay matrix p^-^'^ is obtained 
from p(-^) by replacing q_|_ by — q_ and by K_. The factor K_(. (k_) signifies the top- spin 
analyzing power of the charged lepton. It is equal to one to lowest order in the SM, that is, 
for V — A charged currents. Its value including the order corrections can be extracted 
from the results of [jT^ and turns out to be very close to one: = K = 1 — O-OlStti. 

So far to the building blocks with which the factorizable contributions (which are 
gauge-invariant) to the squared matrix elements of the above reactions are determined at 
NLO. As far as the non-factorizable NLO QCD corrections are concerned which were 



calculated in ref.[|13|], we expect that their contribution to the spin effects is considerably 
smaller than those of the factorizable corrections given below. 

Now we consider the hadronic reactions ([l|) and analyze the following double leptonic 
distribution, 

— — ^ = i(l-Ccos0+cos0_) , (10) 

Of ^^cos0+Jcose_ 4^ +7' 

with Of being the cross section for the channel under consideration. In Eq. ( [T0| ) 0+ (0 ) 
denotes the angle between the direction of flight of the lepton {(.'^) in the t (t) rest 
frame and a reference direction a (b). Different choices will yield different values for 
the coefficient C. Using the general expressions for p, p and the fact that the factorizable 
contributions are of the form Tr[p7?p] we have obtained the following formula for the 
correlation coefficient C in Eq. (p^: 

C = 4K+K_((a-Sf)(b-Sf-)). (11) 
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The expectation value in Eq. dTT] ) is defined with respect to the matrix elements for the 
hadronic production of ttX. It can be expressed in terms of the more familiar double spin 
asymmetry 

zL/ra . vh . ^\ - N(TT)+N(ii)-N(Ti)-N(iT) 

4((a-s,)(b-s,)) - N(TT)+Nai)+N(n)+NaT)' ^^^^ 

where A^(tT) etc. denote the number of pairs with t and t spin parallel - or anti-parallel 
- to a and b, respectively. From Eq. (|T^ one can see that the axes a, b introduced 
in Eq. ( [T0| ) through the angles 0± can be interpreted, as in Eq. to be the spin axes 
of the intermediate tt state within our approximation. This means that the coefficient 
C in Eq. ( [T0| ) reflects spin correlations of the tt intermediate state. Eq. ([11]) holds for 
factorizable contributions to all orders in a^. 

Table 1: Coefficient C of Eq. ([TT]) to leading (LO) and next-to-leading order (NLO) in 



for the spin bases of Eq. (^. The CTEQ parton distribution functions were used and 
we have chosen /jr = /jf = nit = 175 GeV. 





pp at y/s 


= 2TeV 


pp at ^/s 


= 14 TeV 




LO 


NLO 


LO 


NLO 


Chel. 


-0.456 


-0.389 


0.305 


0.311 


Cbeam 


0.910 


0.806 


-0.005 


-0.072 


Coff. 


0.918 


0.813 


-0.027 


-0.089 



Table [I] contains our results for C at leading and next-to-leading order in using the 
parton distribution functions Q] CTEQ5L (LO) and CTEQ5M (NLO). These numbers 
and the results given below were obtained by integrating over the full phase phase. A 
further technical comment is in order here. In order to match with the definition of the 
QCD coupling used in the evolution of the PDF, we have expressed the MS coupling a? 
of 6 flavor QCD by the MS coupling whose evolution is governed by the beta function 
that depends only on the 5 light quark flavors. 

For pp collisions at ^/s = 2 TeV the helicity basis is not the best choice because the 
t, i quarks are only moderately relativistic in this case. Table [T] shows that the dilepton 
spin correlations at the Tevatron are large both in the off-diagonal and in the beam basis. 
These two spin bases yield almost identical results. The QCD corrections decrease the 
LO results for these correlations by about 10%. Since the gg initial state dominates tt 
production with pp collisions at ^/s = 14 TeV the beam and off-diagonal bases are no 
longer useful. Here the helicity basis is a good choice and gives a spin correlation of 
about 30%. In this case the QCD corrections are small. The large difference between the 
LO and NLO results for the correlation in the beam basis at the LHC is due to an almost 
complete cancellation of the contributions from the qq and gg initial state at LO. 

As usual, there are several sources of theoretical errors at fixed order in perturba- 
tion theory. As to the scale uncertainty, the inclusion of the QCD corrections reduces 
the dependence of the ti cross section on the renormalization and factorization scales 
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Table 2: Upper part: Dependence of the correlation coefficients, computed with the PDF 
of ref.[jT^, and /u = /jr = /jf at NLO. Lower part: Correlation coefficients Chei., Cbeamj 
and Coff. at NLO for /jj^ = fjf = nit and different sets of parton distribution functions: 





pp at 


= 2 TeV 


pp aX ^= 14 TeV 




Chel. 


Cbeam 


Coff. 


Chel. 


nit/ 2 


-0.364 


0.774 


0.779 


0.278 


nit 


-0.389 


0.806 


0.813 


0.311 


2mt 


-0.407 


0.829 


0.836 


0.331 


PDF 


Chel. 


Cbeam 


Coff. 


Chel. 


GRV98 


-0.325 


0.734 


0.739 


0.332 


CTEQ5 


-0.389 


0.806 


0.813 


0.311 


MRST98 


-0.417 


0.838 


0.846 


0.315 



significantly. The same is true for the product a^C - see ref.[|T]] for details. To leading 
order in the coefficient C depends only on the factorization scale ^p, while at NLO it 
depends on both scales and ^f. Table ^ shows our NLO results for the three choices 
f^R = I^F = fnt/2,mf, 2mt, again using the PDF of ref.[[T4|]. 

In Table ^ we also compare our results for C obtained with different sets of PDF. 
For pp collisions at ^ = 2 TeV, the spread of the results is larger than the scale un- 
certainty given in the upper part of Table 0. To a considerable extent this is due to an 
important property of C, namely the qq and gg initial states contribute to C with oppo- 
site signs. Therefore the spin correlations are quite sensitive to the relative weights of 
qq and gg initiated tt events. These weights depend in particular on the chosen set of 
PDF. For example, we find the following individual NLO contributions for the helicity, 
beam, and off-diagonal correlation at the upgraded Tevatron: for the GRV98 (MRST98) 
PDF Cgi = -0.443 (-0.486), Cgi = +0.124 (+0.075), C^f^^ = +0.802 (+0.879), 
Cbfam = -0-068 (-0.042), and = +0.810 (+0.889), Cf| = -0.073 (-0.044). This 
suggests that accurate measurements of the dilepton distribution (|TOl), using different spin 
bases, at the upgraded Tevatron could be a useful tool in the effort to improve the knowl- 
edge of the PDF 

The above analysis can be extended to the "lepton-i-jets" and "all jets" decay channels 
in a straightforward fashion [|T7]|. The "lepton-i-jets" channels should be particularly useful 
for detecting tt spin correlations: although one looses top-spin analyzing power one gains 
in statistics and the experimental reconstruction of the t and t rest frames may also be 
facilitated. 

To conclude: we have determined, at NLO in the QCD coupling, the dileptonic angu- 
lar distribution (|TOl) that reflects the degree of correlation between the t and t spins. Our 
results for the Tevatron show that the scale and in particular the PDF uncertainties in the 
prediction of this distribution must be reduced before tt spin correlations can be used in a 
meaningful way to search for relatively small effects of new interactions affecting tt pro- 
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duction and/or decay, but respect, like QCD, parity and CP. For example, the effect of a 
small anomalous chromomagnetic ttg coupling that would lead to deviation from the SM 
predictions at the, say, few percent level would be swamped by these uncertainties. On 
the other hand our Tevatron results should be useful to learn more about the parton distri- 
butions in the proton at high energies. For pp collisions at ^/s =14 TeV the theoretical 
uncertainties in the prediction of this distribution are, fortunately, smaller. This gives rise 
to expectations that top quark spin correlations will play an important role in the precision 
analysis of tt events at the LHC. 
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